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Abstract 

We present our results for inclusive instanton-induced cross-sections in deep-inelastic 
scattering, paying in particular attention to the residual renormalization-scale dependencies. 
A "fiducial" kinematical region in the relevant Bjorken variables is extracted from recent 
lattice simulations of QCD. The integrated instanton-contribution to the cross-section at 
HERA corresponding to this fiducial region is surprisingly large: It is in the 0(100) pb 
range, and thus remarkably close to the recently published experimental upper bounds. 



1. Instantons JT| are non-perturbative gauge field fluctuations. They describe tunnelling tran- 
sitions between degenerate vacua of different topology in non-abelian gauge theories like QCD. 
Correspondingly, instantons and anti-instantons carry an integer topologigal charge Q = 1 and 
Q = —1, respectively, while the usual perturbation theory resides in the sector Q = 0. Unlike the 
latter, instantons induce processes which violate chirality (Q 5 ) in (massless) QCD, in accord 
with the general chiral- anomaly relation. An experimental discovery of instanton-induced events 
would clearly be of basic significance. 

The deep-inelastic regime is distinguished by the fact that here hard instanton-induced processes 
may both be calculated |3]11]] within instanton-perturbation theory and possibly detected experi- 
mentally • As a key feature it has recently been shown || , that in deep-inelastic scattering 
(DIS) the generic hard scale Q cuts off instantons with large size p ^> Q" 1 , over which one has 
no control theoretically. 

In continuation of Ref. [U, where the amplitudes and cross-sections of exclusive partonic sub- 
processes relevant for DIS were calculated, we summarize in the present letter the results of our 
finalized calculations of the various inclusive instanton-induced cross-sections (Sects. 2 and 4). 
A detailed account of our calculations will be published elsewhere 0. The essential new aspect 
as compared to our first estimates is the strong reduction of the residual dependence on the 
renormalization scale resulting from a recalculation based on an improved instanton density []10[ , 
which is renormalization-group invariant at the two-loop level. 

There has been much recent activity in the lattice community to "measure" topological fluctua- 
tions in lattice simulations |TT| of QCD. Being independent of perturbation theory, such simula- 
tions provide "snapshots" of the QCD vacuum including all possible non-perturbative features like 
instantons. They also provide crucial support for important prerequisites of our calculations in 
DIS, like the validity of instanton-perturbation theory and the dilute instanton-gas approximation 
for small instantons of size p < 0(0.3) fm. As a second main point of this letter (Sect. 3), these 
lattice constraints will be exploited and translated into a "fiducial" kinematical region for our 
predictions of the instanton-induced DIS cross-section based on instanton-perturbation theory. 

2. The leading instanton (J)-induced process in the DIS regime of e ± P scattering for large photon 
virtuality Q 2 is illustrated in Fig. [[]. The dashed box emphasizes the so-called instanton- subprocess 
with its own Bjorken variables, 

Q' 2 = -q' 2 >0; x' = ^<l. (1) 

2p ■ q 

As can be inferred from Ref. @ and will be detailed in Ref. ||, the inclusive /-induced cross- 
sectionQ in unpolarized deep-inelastic e ± P scattering can be expressed (in the Bjorken limit) 

as 

m zP ~ V p ' p o- {I) (r' O' 2 ) (2) 

dx> dQ' 2 ~ fr; dx< dQ' 2 p,p[ ' ^ 1 { ) 

where p' = q',q' denotes the virtual (anti-)quarks entering the J-subprocess from the photon side 
and p = q,q,g denotes the target partons (c.f. Fig. [I]). The differential luminosity dC^) p , account- 
ing for the number of p'p collisions per eP collision, has a convolution-like structure || , involving 
{^Bjj 2/Bj, x}-integrations over the target-parton density, f p (xBj/x, . . .), the 7*-flux, P 7 *(?/Bj), and 

1 A sum over /-induced (AQ5 = 2n/) and anti-instanton (J)-induced (AQ5 = —2nf) processes is always implied 
by the superscript (/) at cross-sections. 
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the known 



flux P { J\x/x', 



of the parton p' in the /-background. We shall display the 



explicit form of the differential luminosity in Sect. 4. 
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Figure 1: The leading instanton- induced process in the DIS regime of e ± P scattering. 

The simple relation (||) between da e p and cx^, derived within /-perturbation theory |7|,|| in the 
Bjorken limit, is actually much less obvious than an inspection of the grossly oversimplified Fig. [I] 
may suggest. The derivation proceeds in two steps: i) Using the Feynman rules of /-perturbation 
theory in momentum space, one calculates the (manifestly gauge invariant) inclusive eP cross- 
section, with the current quark being not a free parton, but rather described by the (complicated) 
quark propagator in the /-background as in Ref. ii) One independently writes down the 
(gauge invariant) expression for the total cross-section a^ p with an off-shell external parton p'. 
In the Bjorken limit, when certain non-planar contributions may be neglected, one then finds by 
comparison of i) and ii) the form (0) of the inclusive /-induced eP cross-section along with the 
explicit expression for the flux-factor Pfp . 

While a more detailed description of this rather involved calculation has to be deferred else- 
where 101, let us summarize next the state of the art evaluation of the /-subprocess total cross- 
Q' 2 ), which contains most of the crucial instanton-dynamics. 

We start with the /-subprocess total cross-sections (here only for the dominating case of a target 
gluon) in a form |7]|| still exhibiting the complicated integrations over collective coordinates^ 
(/-sizes p,p, ...), 
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3a s (p r ) {{p + q'Yf 2 



(3) 



where p' = q', q'. The most important quantities entering Eq. @ are: 

2 For brevity, we display the cross-section (||) already after the (saddle-point) integration over the /-colour 
orientations. The function O, whose explicit form will be specified below, accounts for that. 
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The I-density D(p,p r ), which has the general form 

D{p ,^ d ( J* V«pf-^|— (4) 




with p r denoting the renormalization scale and [|l(J 

MltJ* A 2S 12,3„^, (5) 
in terms of the QCD /^-function coefficients, /3 = 11 — /3i = 102 — |n/. The power 



A)Ai — A2 makes the /-density renormalizat ion-group invariant at the two-loop level |[L0| , 
(l/D)dD/dp r = 0(a 2 s ), in contrast to the original one-loop expression 0, corresponding 
to Ai = 1 and A2 = 0, with (1 / D)dD / dp r = 0(a s ). The constant d is scheme-dependent; 
in the MS-scheme it is given by || d = dexp[-3C 2 + n f C 3 )/2, with C x = 0.46628, 
C 2 = 1.51137, and C 3 = 0.29175. 

Note that the large, positive power of p in the /-density (Q) would make the integrations 
over the /-sizes in Eq. (||) infrared divergent without 

• the form factors Ki(Q'p(jj)): For large Q'p(p), the virtuality Q' of the internal quark p' in 
Fig. [1] provides an exponential cut-off, K\{Q' p) oc exp(— Q'p); m the integrations over the 
/-sizes Q. These form factors where shown to arise naturally in step i) above, which is 
manifestly gauge invariant ||. In step ii) one has to adopt a gauge- invariant definition of 
the p'p cross-section, since the incoming parton p' is off-shell |[L4|| . Then one obtains exactly 
the Bessel-K form factors [Bj, unlike naive, not manifestly gauge-invariant definitions which 
lead in addition to these well-defined contributions to unphysical ones suffering from infrared 
divergent /-size integrations |T5| , P^| ]. 

• The functions Q and u (along with the integration over R^) summarize the effects of final- 
state gluons (Q) and final-state quarks (uj). The function Q, appearing in the exponent with 
a large numerical coefficient, 4tt/q: s , and u, occuring with a high power, 2nj — 1, call for a 
precise evaluation. Hence, let us turn next to describing their state-of-the-art evaluation. 

It is very instructive to consider two alternative interpretations of the functions Q, u, and the 
integration variable R^. 

• Total cross-section via summation of exclusive cross-sections: 

This is the cleanest and most straightforward method to arrive at the total cross-section. 
In this case one starts with the familiar representation of the (^•'-function associated with 
energy-momentum conservation, 

(2 tt) 4 5 {i \p + q' - £ h) = J d A R exp [i (p + q' - £ h) ■ R]. (6) 

i i 

The phase-space integration over the final-state gluons/quarks is then performed by means 
of the basic formula 

/^^.PH^.S]^-^, (7) 
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with the help of which one finds 
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-R 2 + ieR 
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3/2 



+ . . . . 



(9) 



The interpretation of the various terms contributing to the perturbative expansion of Q 
in Eq. (|S[) is illustrated in Fig. ^| (left): The first term takes into account the summation 
and exponentiation of the leading-order gluon emission, whereas the second term originates 
from the summation and exponentiation of interference terms between the leading-order 
gluon emission and the gluon-propagator correction. The first term contributing to the 
perturbative expansion of uj in Eq. @ just corresponds to the leading-order quark emission. 
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Figure 2: Left: Graphical interpretation of an exclusive /-induced q'g cross-section in terms of 
perturbation theory in the /-background. (Curly) lines ending at blobs denote LSZ-amputated 
quark zero modes (classical /-gauge fields). Curly lines connected by an ellipse denote LSZ- 
amputated gauge-field propagators in the /-background. Right: The //-valley action correspond- 
ing to the most attractive //-colour orientation. 



In summary, the perturbative approach based on the exclusive amplitudes, as calculated 
within /-perturbation theory, yields the essential functions fl and uj as asymptotic expan- 
sions for small pp/R 2 . Since p and p are conjugate to the virtuality Q' and R is conjugate 
to the total momentum of the /-subprocess, p + q', (c.f. Eq. @ and Fig. [I]), we expect 
qualitatively 

p ~ p ~ l/Q' and R 2 ~ l/(p + q'f pp/R 2 ~ (p + q'f/Q' 2 = 1/x' - 1. (10) 

Thus, strict /-perturbation theory for the total cross-section is only applicable for not too 
small x' . 

Total cross-section via optical theorem and 1 1 -valley method: 



In this approach, one evaluates |18|] the total cross-section from the imaginary part of the 



forward elastic scattering amplitude induced by the instanton-anti-instanton (//)-valley 
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background, A^ T \ In this case, stands for the separation between / and /, and VI is 
identified with the interaction between I and /, 

n~ s^co-i. (n) 

In the valley approximation, the //-valley action, S^ 11 ^ = S[A^^], is restricted by con- 
formal invariance to depend only on the "conformal separation" Jl9| 



-fl» + kfi. +g + g (12 ) 

PP P P 

and its functional form is explicitly known [|15|,f2(]] (Fig- @ (right)). 

Note that for all separations £, the interaction between / and / is attractive (c.f. Fig. |2] 
(right)): The //-valley corresponds to a configuration of steepest descent interpolating 
between an infinitely separated /// pair and a strongly overlapping one, annihilating to the 
perturbative vacuum. 

Analogously, the function u is now identified with the fermionic overlap integral for which 



an integral representation was found in Ref. pi], which we were able to perform analytically, 



c(0 = 2 /i ( 5 *;4; 1 - \) ; z = I (e + J? - 4) . (13) 



3/2 ' * y 2 ' 2' ' W ~ 2 

Finally, the function Cl arising from the integration over the relative //-colour orientations 
has been estimated in Ref. [[0| by assuming for simplicity an orientation dependence of the 



valley action^ corresponding to a dipole-dipole interaction |22 

fi<€)*e^. d4) 

The leading terms in the asymptotic expansions of the //-interaction ( ITT] ) and the fermionic 
overlap ( Ti3|) for large conformal separation, 

m = ~| + ^(ln(0/a), w(0 = + 0(ln(O/^ 7/2 ), (15) 

exactly reproduce the known perturbative results (§) and (|]) for small pp/R 2 . This il- 
lustrates the power of the //-valley method to effectively sum up the gluonic final-state 
tree-graph corrections to the leading semi-classical result^ |T7| . 



We shall thus take the valley expressions for Q, u and Cl, Eqs. (|TT]), (0) and (0), to smoothly 
extrapolate somewhat beyond strict /-perturbation theory. 

Let us add, however, that the full content of the valley approximation is not essential in this 
context. It is mostly the shift in the expansion variable 

PP _ PP -1 (16) 



-R 2 + ieflp -R 2 + ieR + p 2 + p 2 £ 



3 The saddle-point corresponds to the most-attractive //-orientation. We have checked |j| that taking into 
account the exact orientation dependence of the valley action [ pp| gives numerically a very similar result. 

4 Some initial-state and initial-state - final-state corrections might exponentiate as well pg| ]. These are not taken 
into account by the valley action. 
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making the leading terms (|15D qualitatively adequate down to fairly small £( >3), in contrast to 
the strict i-perturbative expansions (^) and (0). 

The collective coordinate integration in the cross-section (|3|) is perfectly suited for a saddle-point 
evaluation. To this end, we collect the most-relevant factors in Eq. (|3|) in the following effective 
exponent, 

-r = i(p + <t)'R- Q'(p + p)- (^^-A^^^ . (17) 



In arriving at Eq. ( 17 ) we have used the asymptotic form Ki(Q'p) oc exp [—Q'p] for the Bessel-K 
functions, anticipating that, for small a s (p r ), the dominant contribution to Eq. (|3|) will come 
from the region Q' p(p) ^> 1. Note that the parameters Ai and A2 allow us to trace the impact 
of the two-loop improvement of the /-density, with the one-loop expression || corresponding to 
Ai = 1 and A 2 = 0. This is to be contrasted with previous related studies which either ignored 
the crucial renormalization-scale dependences altogether [|T5l , |T6l or were still too crude for a study 
of the associated uncertainties [0] (c.f. also Fig. |3] (left) below). 

The corresponding saddle-point in R^, p and p is most easily found in the p'g centre-of-mass 
(cm.) system. One finds R* = (—ip*y/£* — 2, 0) and p* = p*, where £* and p* are the solutions 
of the following saddle-point equations, 



Q'P* - ( — ~7~T ~ 2 A i A) In ( P *Pr)) ^l!^ = 0, (1* 



1 1-x' 



2 V x' 



e-2-1 \Q l p* + A 1 p S {II \C)-^2 = 0. (19) 



Upon evaluating the integrand in Eq. (|3|) at the saddle-point and taking into account the in- 
tegration over the (Gaussian) fluctuations about the saddle-pointf], we may finally express the 
cross-section entirely in terms of v* = Q'p* and £*, 

Q' 2 4l = d 2 ^^((C + 2K 2 + 45(5 - 2v*)) ~ 2) V2 ^(D 2 "^ 1 (20) 



(e - 2f v 

x = — 



3„,* 5 



(v* - 5) 9 / 2 + 2)v* -4Sj±(S-v*- 2D(S)) 2 + 5(5 - ( In [^== 



47T 



19/2 



exp 



4 " -^(D-2fl-lnf^U' 



H'2„W - 32 q 8 (/x r ) 1 yM 1 ^ (J) (J) _ n r w (/) / 91 x 

where we have introduced the shorthands 

5(f) = A 1 (3 S^(C) ~ A 2 , D(f(C)) = d ^ /(D. (22) 

5 We have checked that our result for the Gaussian integrations coincides, for the one-loop case (Ai = 1, A2 = 0), 
with the corresponding result quoted in Ref. [p^| . 
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For completeness, we have listed also the corresponding expression for the q'q cross-section (|2l|) . 

What remains is to solve the saddle-point equations, (|T8| ) and (|19|). An analytical solution || 
in the asymptotic regime a s (Q') — > 0, x' and fi r /Q' fixed, confirms our qualitative expectations 
(|T(]|). In particular one finds, asymptotically, R*/p* = — 2 — > 2/^/l/V — 1. However, for 
experimentally accessible values of the virtuality Q', the corrections to the asymptotic result are 
quite large and the corresponding analytical expressions complicated. Hence, we only present 
here the results corresponding to a numerical solution of the saddle-point equations. 
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Figure 3: Renormalization-scale dependences of the J-subprocess cross-sections for a target 
gluon, Eq. (pOD, (left) and a target quark, Eq. (|2~lD, both for rij = 3. 

In Fig. ^, we display the residual renormalization-scale dependencies of the J-subprocess cross- 
sections (|20|) and (|21|) over a /arg'e range of \x T jQ' . Apparently, we have achieved great progress 



in stability and hence predictivity by using the two-loop renormalization-group invariant form of 
the /-density D(p,fi r ) from Eqs. (|4]) and ([5]): The residual dependence on the renormalization 
scale \x r turns out to be strongly reduced as compared to the one-loop case (Ai = 1, A 2 = 0). 

Intuitivelely one may expect [^|,|TB],^| \i r ~ l/(p) ~ Q'/ Po — O(0.1) Q'. Indeed, this guess turns out 
to match quite well our actual choice of the "best" scale, /i r = 0.15 Q', for which do^) g jdp, r ~ 
(c.f. Fig. |] (left)). We also note that the cross-sections for a target gluon (Fig. [| (left)) are 
about two orders of magnitudes larger than the cross-sections for a target quark (Fig. |3] (right)). 
Henceforth, the latter are neglected. 

Our quantitative results on ajp g are shown in detail in Fig. f|, both as functions of Q' 2 (left) and of 
x' (right). The dotted curves indicating lines of constant p* (left) and of constant R*/p* (right) 
nicely illustrate the qualitative relations (|T0|): For growing Q' 2 and fixed x', smaller and smaller 
instantons, p* ~ l/Q', are probed and the cross-sections decrease rapidly, mainly because of the 
large powers of p in the /-density (f|). For decreasing x' and fixed Q' 2 , on the other hand, the II- 
separation R* in units of the /-size p* decreases and the cross-section increases dramatically. In 
the language of the //-valley method the latter originates mainly from the attractive interaction 
between instantons and anti-instantons. 

3. We have seen that the collective coordinate integrals in (|3|) are dominated by a single, calculable 
saddle-point (p*,R*/p*), in one-to-one relation to the conjugate momentum variables (Q',x'). 
This effective one-to-one mapping of the conjugate /-variables allows for the following important 
strategy: We may determine quantitatively the range of validity of /-perturbation theory and the 
dilute /-gas approximation in the instanton collective coordinates (p < p ma ,x,R/p > (R/p)mm) 
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Figure 4: The instanton-subprocess cross-section (p0| ), for rif = 3, both as functions of Q' 2 
(left) and x' (right). The dotted curves are lines of constant J-size p* (left) and of constant 
//-separation R* in units of the /-size p* (right). 



from recent (non-perturbative) lattice simulations of QCD and translate the resulting constraints 
via the mentioned one-to-one relations into a "fiducial" kinematical region (Q' > Q' min , x' > x' min ). 

In lattice simulations 4d-Euclidean space-time is made discrete; specifically, the "data" from the 
UKQCD collaboration [p4|, which we shall use here, involve a lattice spacing a = 0.055 — 0.1 
fm and a volume V = i S p ace • /time = [16 3 ■ 48 — 32 3 ■ 64] a 4 . In principle, such a lattice allows to 
study the properties of an ensemble of (anti-)instantons with sizes a < p < V 1 ^. However, in 
order to make instanton effects visible, a certain "cooling" procedure has to be applied first. It 
is designed to filter out (dominating) fluctuations of short wavelength 0(a), while affecting the 
topological fluctuations of much longer wavelength p ^> a comparatively little. For a discussion of 
lattice-specific caveats, like possible lattice artefacts and the dependence of results on "cooling" 
etc., see Refs. [0.^4] . 
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Figure 5: Support for the validity of /-perturbation theory for the /-density D(p) (left) and the 
dilute /-gas approximation (right) for p < p max ~ 0.3 fm from recent lattice data [j24| . 

The first important quantity of interest, entering /-induced cross-sections (c.f. Eq. (^), is the 
I-density D(p), Eq. (f|). This power law, D(p)\ n =0 oc p 6 , of /-perturbation theory is confronted 
in Fig. [I (left) with recent lattice "data", which strongly suggests semi-classical /-perturbation 
theory to be valid for p <p max — 0.3 fm. Next, consider the square of the total topological charge, 
Q 2 = (n — n) 2 , along with the total number of charges, N tot = n + n. For a dilute gas, the number 
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fluctuations are poissonian and correlations among the n and n distributions absent, implying 
(Q 2 /N tot ) = 1. From Fig. ^ (right), it is apparent that this relation, characterizing the validity of 
the dilute /-gas approximation, is well satisfied for sufficiently small instantons. Again, we find 
p max ~ 0.3 fm, quite independent of the number of cooling sweeps. For increasing p max ^0.3 fm, 
the ratio (Q 2 /N tot ) rapidly and strongly deviates from one. 

Crucial information about a second quantity of interest, the 1 1 -interaction, may be obtained as 
well from the lattice [0,E3]. Quite generally, it is found that the semi-classical attraction for 



large R 2 / (pp) turns into a non-perturbative repulsion for smaller separations in units of the sizes, 
such that in vacuum^] (R 2 /(pp)) = 0(1). Thus it seems a reasonable extrapolation to use the 
attractive, semi-classical valley result for the J J- interaction Q, Eq. ([11]), down to a minimum 
conformal separation £ min ~ 3, corresponding to (-R*/P*)mm — 1- 

Finally, by means of the discussed saddle-point translation, these lattice constraints may be turned 
into a "fiducial" kinematical region for our cross-section predictions in DIS (c.f. Fig. |]), 



P* < P^ ax ^ 0.3 fm; ) { Q' > Q' ~ 8 GeV 



P — \P 




£ > M£) - I I I x' > xL„ ^ 0.35 

min 



(23) 



Unlike DIS, where only small instantons are probed, in the /-liquid model of Ref. |27j more 



emphasis is placed on the physics associated with larger instantons. For /-ensembles includ- 



ing also larger /-sizes >0.3 fm, the various recent lattice results [|TTJ[y.[2^.f2(| do not, however, 
unanimously support the liquid picture. 

4. Experimentally, in deep inelastic eP scattering at HERA, the cuts (|23"D must be implemented 
via a (Q',x') reconstruction from the final-state momenta and topology ||, while theoretically, 
they are incorporated into our /-event generator |J "QCDINS 1.6.0" and the resulting prediction 
of the /-induced cross-section in DIS at HERA. The latter is connected to the /-subprocess 
cross-sections a^, by the differential p'p luminosity || (c.f. Eq. (Q)), 

We| 'j «fc J d^'p-M PA w)^,...)/^,...). (24) 
ax' dQ' 2 b x' J x J x B j J yBj x % 

■^Bj min *^Bj min ^Bj min 

Here S (~ 9 • 10 4 GeV 2 for HERA) denotes the cm. energy squared of the eP collision, e 2 , is 
the electric charge squared of the current (anti-) quark in units of the electric charge squared, 
e 2 = Ana, and P 7 » denotes the familiar Weizsacker- Williams-type photon flux, 

P 7 ,(y B j) = (l + (l-2/Bj) 2 )A/Bj, (25) 

with |/Bj = Q 2 j \Sx-Qyj- Furthermore, f p (x^/x, . . .) denotes the density of the target parton p in 
the proton, with the dots standing for the factorization scale, and, finally, the factor pjy accounts 
for the flux of virtual (anti-) quarks p' in the /-background entering the /-induced p'p-subprocess 
from the photon side |7||| (c.f. Fig. [j]), 




16 7T 3 x' \ x x' Q 2 



Published ratios range from (R)/(p) ~ 0.83 @, (R/{p + p)} ~ 0.59 |§ to (R/(p + p)} ~ 1 || 
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The /-induced cross-section in DIS at HERA, c>h ERA , subject to kinematical cuts (scbj > ^Bjmin; 
2/Bjmax > Itej > Z/Bjmin; > a^; Q' > Q' min ) , is then obtained by integrating Eq. (g) over the 
appropriate range of x' and Q' 2 . 

Let us point out that the factorization of the eP cross-section for fixed x' and Q' 2 into a sum 
of differential luminosities and i-subprocess cross-sections, is essential for the possibility to place 
different cuts on the Bjorken variables of the eP and the p'p system, respectively Of particular 
interest is XBjmin <C x' min . Such cuts permit to explore essentially the full accessible (xbj,x, Q 2 ) 
range in DIS at HERA, down to (10~ 3 , 10~ 3 , 10 GeV 2 ), say. By placing in this region the additional 
cuts dffi) on (x', Q' 2 ), /-searches benefit from the high statistics at small x-Qj, while the theoretical 
control is retained over the /-dynamics. 

In Ref. [||], on the other hand, (only) the infrared safe pieces of the /-induced contributions to 
the parton-structure functions, J-2p(x, Q 2 ), were estimated in one step by means of configuration 
space techniques. In our momentum space language, the authors have implicitly integrated over 
Q' 2 and x', with xej < x < x' < 1. Hence, the results of Ref. || can only be applied |5|] to 
relatively large XBjmin = ^min ~ 0.35. 

Nevertheless, as an important check of our calculations, we have also calculated the infrared-safe, 
/-induced contributions to the parton-structure functions, by integrating our asymptotic results 
over x' and Q' 2 and retaining only the contributions from the upper Q /2 -integration limit cx Q 2 . 
Within the common range of validity |J of various employed approximations, we find perfect 
agreement with the gluon-structure function quoted in Ref. @. 




Figure 6: Instanton- induced cross-section at HERA [nj = 3). 

Fig. H displays the finalized /-induced cross-section at HERA, as function of the cuts x' min and 
Qmi n , as obtained with the new release "QCDINS 1.6.0" of our /-event generator. Only the target 
gluon contribution has been taken into account. For the minimal cuts (^) extracted from lattice 
simulations, we specifically obtain 

^HERA<y > 0.35, Q' > 8 GeV) ~ 126 pb; for x B] > 10~ 3 ; 0.9 > y B] > 0.1. (27) 

Hence, with the total luminosity accumulated by experiments at HERA, C = O(80) pb -1 , there 
should be already (9(10 4 ) /-induced events from this kinematical region on tape. Note also that 



11 



the cross-section quoted in Eq. ( p7|) corresponds to a fraction of /-induced to normal DIS (nDIS) 
events of 

f {1) = -J^k = °M % ! for ^bj > 10- 3 ; 0.9 > y Bj > 0.1. (28) 
°"hera 



This is remarkably close to the published upper limits on the fraction of /-induced events [ f29f , 
which are also on the one percent level. 

There are still a number of significant uncertainties in our cross-section estimate. For fixed Q'- 
and x'-cuts, one of the dominant uncertainties arises from the experimental uncertainty in the 
QCD scale A. We used in the two-loop expression for a s with rif = 3 massless flavours the value 
= 282 MeV, corresponding to the central value of the DIS average for n f = 4, Ag, = 234 

MeV [p5| . If we change A^± within the allowed range, ~ ±65 MeV, the cross-section (p7|) varies 
between 26 pb and 426 pb. Minor uncertainties are associated with the residual renormalization- 
scale dependence (c.f. Fig. |3|) and the choice of the factorization scale. Upon varying the latter 
by an order of magnitude, the changes are in the (9(20) % range. 

By far the most dominant uncertainty arises, however, from the unknown boundaries of the 
fiducial region in (x', Q') (c.f. Fig. |6|). Here, the constraints from lattice simulations are extremely 
valuable for making concrete predictions. 
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